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Abstract 

Recently, spin-selective radical pair reactions have been studied using concepts from quantum 
measurement theory. In this Article, we show that the approach taken by Kominis (Physical Review 



O ■ 

. E, 83, 2011, 056118) leads to erroneous results due to a problematic treatment of quantum jumps. 

Correct consideration of quantum jumps leads to either the traditional master equation or the 
Jones-Hore master equation. 
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I. INTRODUCTION 



Spin-selective radical-ion pair reactions are at the core of spin chemistry. The phe- 
nomenological master equation [1] (1) is often used to model spin-selective reactions. 

d { t =^[H,p]- k f{Qs,p}- k f{Q T , P } (i) 

where p is the density matrix of the reactants, {} denotes the anti-commutator, H is 
the Hamiltonian describing unitary evolution of the radical pair, k$ and kr are the rates of 
reaction through singlet and triplet channels respectively, and Qs and Qt are singlet and 
triplet projection operators respectively. Recently, based on quantum measurement theory, 
Kominis J2, ^ derived equation (2) and quantum jump equations (3) and (4) 

^ = ~i[H, P] - ^^(pQs + QsP ~ 2QspQs) (2) 

p nr describes the state of the radical pairs before recombination. Equation (2) is derived 
by making the analogy between spin-selective chemical reaction and electron tunneling in 
a quantum dot. The quantum non-demolition measurement of the quantum dot causes the 
wavefunction to collapse periodically, thus pure decoherence without energy dissipation is 
observed. Kominis then introduced quantum jumps (3) and (4), where ps is the probability 
of reacting down the singlet channel and Pt is the probability of reacting down the triplet 
channel, to represent chemical reaction. 

Ps = k s (Qs) dt (3) 

Pt = k T (Q T ) dt (4) 

This an erroneous application of the Bohr-Einstein quantum jump |4| and leads to misin- 
terpretation of the "no jump" event [5| which will be analysed below. Subsequently, Kominis 



6] 



introduced a revised master equation based on a phenomenological interpolation between 
the "maximal coherence case" and the "minimal coherence case" without correcting the un- 
derlying physics in (3) and (4). Rephrasing the argument presented, singlet-triplet coherence 
is measured by the parameter p co h- 

h = Tr{p ST PTs} 
coh Tr{p ss }Tr{p TT } 



2 



where p ST = QspQt, Pts = QtpQs, Pss = QspQs and p TT = QtpQt • The evolution of 
the density matrix is given by 



dp tit i k s + k T ,dp incoh dp coh 

— = -i[H, p\ {pQs + Qsp - ZQspQs) - (1 - pcoh)—j- t pcoh ~fo~ ( 6 ) 



dpincoh = k s dtQspQs + k T dtQ T pQ T (7) 

dpcoh = {ksTr{Q s p} + k T Tr{Q T p])dt - (8) 

Tr{p\ 

The interpolation parameter p co h is introduced "by hand" and the limiting cases of the 
resulting equation will be commented on below. This Article will also show how consistent 
derivation of the master equation points towards (9), the Jones-Hore master equation Q or 
the phenomenological master equation. 

f = -i[H, P \ - 'fiQs, P } - ^{Q T , P } - ^^(pQs + QsP - 2Q sP Q s ) (9) 

A microscopic derivation of (1) has been reported by Ivanov et al. [8j. Appendix A 
contains the full microscopic derivation of (9). 



II. COMMENT ON THE KOMINIS MASTER EQUATION 
A. Problem with "no jump" events 

The quantum jump equations (3) and (4) fail to capture the physical significance of a 
no-jump event. In a qualitative sense, a no-jump event represents either the wavefunction 
is still in {S, T} subspace or, importantly, the fact that the radical pair has already reacted 

Hi- If one starts observing the system at time t and does not see the system executing 
a quantum jump after a very long time, one should conclude that the system has already 
jumped at some prior time before the observation rather then believing that the system is 
yet to jump. To capture this, a non-Hermitian term must be added to the Hamiltonian 

H eff = H-i^.\S)(S\-i^\T)(T\ (10) 

Putting into (2) the non-Hermitian Hamiltonian (10), which corresponds to quantum 
jump equations (3) and (4), and only considering the S-T subspace 



dp dp reac n on dp n 



dt dt dt 

= -i[H -i^\S)(S\-i J ^\T){T\,p\- ^±^{ P Q S + Q sP - 2Q sP Qs) 

= -i[H, p] - ^{Q Sl p} - ^{Q T , p} - ^±^{pQ s + Q sP - 2Q sP Qs) 

and the Jones-Hore master equation [7] is recovered. 

The role of the non-Hermitian Hamiltonian can be seen more clearly when one writes the 
jump-free evolution of the system explicitly. When ||ff e j/|| dt <C 1 

m + dt)) = e- iH °ff dt \i>{t)} 

« (i - iH eff dt) \m) (ii) 

ignoring second order terms in dt 



\tfj(t + dt)\ 2 = (lP(t)\(l+iH+ ff dt)(l- iH eff dt) \lP(t)) 

= 1-Sp (12) 

where dp is the decrease in the norm of the wavefunction which is compensated by the jump 
to products. Identifying dp as the probability of a quantum jump 

5 P = idt(m\H e} f-Hj ff \m) 
= dt(k s m)\ Qs \m) + k T m)\ q t \m)) 

= ksdt(Q s ) + k T dt(Q T ) (13) 

This shows qualitatively that although the jump operators "fill up" the states representing 
the reaction products (Sq and T ), the populations of the radical pair S and T states are 
only "removed" by the non-Hermitian Hamiltonian. 

From an algebraic perspective, any completely positive Markovian evolution can be writ- 
ten in the Lindblad form 

% =-i{H, A+ D( P ) 

= -i[ff, /3 ]+X>pL+-i{L+L,,p} (14) 

i 
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where Li are Lindblad operators. The formal solution to (14) can be written as 

Pit) = e Lt P (Q) (15) 
L is a superoperator that can be written trivially in terms of the jump superoperator J 

L = J + L-J (16) 
using an identity for superoperators {ll| 

^ / dx k dx fc _i... / dx l e h{x - Xk) be h{Xk - Xk -^b...be ixi (17) 



OO 

,(a+6)a; 



fc=0 

e ()+l-i) tp(0) = y> /"'^ f tk dt^... r^e^-^^'^ie^-^^-^-^i.Je^-^V^) 
fc=0 ^ ^ ^ 

(18) 

Reading (18) from right to left, L — J can be interpreted as the "between jump" evolution 
of the system and at t\ the system experiences a first jump, followed by a period of between 
jump evolution and so on. J, the jump operator, can be identified with terms bilinear in p. 

J=j>.l+ ( 19 ) 

i 

L = -iH* + *iH + Y j ^L 1 * L+ - h+L t m~m L+L^j (20) 
L-J = - iH m + miH-J2 (l L t L i # +\ # L i L i) 

= -i(H • (^Ltu • +1 • Ltu\ ) 

i ^ ' 

= -AH- l -Y L t L ^^ ( 21 ) 

i 

Identifying the fact that if a quantum jump is not executed between time interval t and 
t + dt, e ( - L ~ J ^ dt propagates the state p(t) — » p(t + dt). 

\^{t + dt)) = e - iH ^f dt \i/}(t)) (22) 

where 

H eff = H- l -Y,LiL l (23) 
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Substituting the jump operators corresponding to (3) and (4), which will be discussed 
below, into (23) we obtained (10). This shows clearly how the non-Hermitian Hamiltonian 
is central to the quantum jump formalism. It is Kominis's neglect of this quantity that led 
to erroneous results. 

B. Expression for dpi nco h 

If there are no coherences, p = QspQs+QtpQt and equation (7) follows trivially. Writing 
(6) in the limit p coh = 



The Jones-Hore equation (9) also reduces to same equation. 
C. Expression for dp co h 

The expression for dp co h suggests that the total density matrix is removed at the combined 
rate at which the reactants are transformed to product. This approach is no longer state 
selective as we know that \S) —> \T ) and \T) — > \S ) transitions are forbidden. Hence a 
differential rate between singlet recombination and triplet recombination should manifest 
itself in different rates at which the singlet and triplet states are depopulated. This suggests 
that the expression for dp co h is correct only if ks = kr- In the limit p co h = 1 and ks = kx = k, 
(6) reads 



dp 
~dt 



i[H,p] 




t[H,p] 



dp 
~dt 



i[H, p] - k(pQ s + Q sP - 2Q sP Qs) - (kTr{Q sP } + kTr{Q T p}) 



Tr{p} 



P 



(24) 



Knowing that 



Tr{p} = Tr{Q sP } + Tr{Q T p} 



(25) 



we find 



-£ = -i[H, p] - k(pQ s + Q s p - 2Q sP Qs) - kp 
(26) is exactly equal to the Jones-Hore equation (9) in the same limit. 
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(26) 



D. Unphysical Prediction of the Kominis Master Equation 



The Kominis master equation predicts that in the absence of singlet-triplet interconversion 
and kj- = 0, starting from a totally coherent mixture of S and T one will have triplet 
population equal to 0.25 after reaction J(5|. 

Kominis attempted to defend the unexpected fall in (Qt) by analysing a single molecule 
trajectory. Paraphrasing his argument, at t — 0, the radical pair can either react with 
probability p r = ksdt (Qs) — or n °t react with probability p nr = 1 — Conditioned 
on the fact that the radical pairs do not react, measurement at rate ^ causes the singlet 



projection to occur with probability qs = ks ^ s ' )dt = , triplet projection to occur with 



2 

4 " , triplet pro; 

with probability q? = ks ^^ dt = hs^t an( j probability of no projection qo = 1 — As 
pure singlet will react eventually and pure triplet will never react, summing up pure singlet 
produced by the measurement and singlet product will give the total singlet yield. 

Ys = (Pr + PnrQs) + PurQaiPr + PnrQs) + (Pnr?fl) 2 (Pr + PnrQs) ■ ■ ■ (27) 



ksdt 
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Pnrqo = 1 - « 1 - k s dt (28) 



ksdt ( k s dt\ k s 3k s dt 
Pr +p nr qs = — + ( 1 2~ J T ~ 4~~ ^ ' 



By summing I5, Kominis obtains 



n=0 

The error in this summing procedure is the fact that the whole density matrix is removed 
with probability p r = ksdt (Qs) to form the singlet product in the first step. This leads to 
non-conservation of spin angular momentum. 

Kominis attempted to circumvent this lack of conservation of spin angular momentum 
by invoking a time averaged p co h- He argues that 

{(Tr{ps T (t)pTs(t + r)})) 

Pcoh{t - ^—f 777T7 1 W L 

Tr{p ss }Tr{p TT } 

PTs(t + r) = e- lHT p TS (t)e iHT (32) 

where ((. . .)) indicates time average over r, with r being larger than the inverse S-T 
energy difference and smaller than the characteristic time-scale of the reaction. Assuming 
the S-T energy separation is J, as e~ lJt rotates rapidly around the complex plane 



Pcoh{t - —- r— — ; PS U (66 

Tr{p ss }Tr{p TT } 

However, (31) is physically questionable. There is no necessary relationship between the 
S-T energy difference and the characteristic time-scale of the reaction, both are system- 
dependent parameters. 

Furthermore, the validity of (32) is dubious. By integrating the Liouville-von Neuman 
equation 

p{t + r) = e~ iHl p(t)e lHT (34) 



p TS (t + r) = Q T e- %UT p(t)e tMT Q s (35) 

only when [H, Q$] = and [H, Qt] = 0, i.e. in the absence of singlet-triplet interconversion, 
can the order of the exponentials be reversed, so that (32) is recovered. 



E. Radical-ion-pair reactions and the optical double slit 

The introduction of p co h and the analogy with the optical double slit experiment suggests 
that there is singlet interference to triplet product and triplet interference to singlet product. 
In an optical double slit, the field incident on the screen at position at time t is the 
superposition of the fields from the two slits 

E M (!>,t) = El +) (l>,t) + E^ ) (!>,t) (36) 

This is because the two slits are physically identical in an optical double slit experiment 
and a single photon can be diffracted by both slits. In a chemical reaction, spin angular 
momentum conservation demands that only singlet reactant can enter the "singlet reaction 
slit" and only triplet reactant can enter the "triplet reaction slit". A superposition state 
cannot react, the coherence is destroyed as the starting state for the jump is selected from 
among the stationary states represented in the superposition. 
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III. CONSISTENT DERIVATION OF MASTER EQUATION USING QUANTUM 
JUMP APPROACHES 



A. Consistent Derivation of Jones-Hore Master Equation 

Treating quantum measurement and chemical reaction as quantum jump processes gives 
a physical interpretation of the Jones-Hore master equation. Using jump operators (37) - 
(40) , the master equation (9) can be recovered. 

(37) 

(38) 
(39) 
(40) 

(37) and (38) are the same as (3) and (4) and correspond to quantum jumps to singlet 
and triplet product respectively. (39) and (40) correspond to a strong measurement of the 
system, the same intuition that lead to equation (2). The evolution of the wavefunction, 
\ip(t)) ->■ \ip(t + dt)) is described by 

\^,{t + dt)) = |yv ;/ p = 1 - k s dt - k T dt 



Jl = 


\S ){S 


J 2 = 


\To)(T 


■h = 


\T)(T\ 


Ja = 


\S)(S\ 





:,H+ f dt e 


\i(j(t + dt)) - 


= \S ) ■ 


\il>(t + dt)) 


= |T> = 


\ll>(t+dt)) = 


= |T > ^ 


\i/>(t + dt)) 





■h \m) 

h W)) 



P = k s (Qs) dt 
P = k s (Qt) dt 
p = k T (Q T ) dt 

p = k T (Q s ) dt (41) 



The form of the non-Hermitian effective Hamiltonian describing both quantum measure 
ment and chemical reaction (42) follows from (23) 

i 



H eff = H ~ 7} XI L i Li 



= H - % - (k s Q s + k T Q s + k s Q T + k T Q T ) 

= H-r k -^L (42) 



The interpretation of (41) is that with probability ksdt, the wavefunction attempts to 
react via the singlet channel and with probability kxdt, the wavefunction attempts to react 
via the triplet channel. The singlet channel measures the wavefunction and with proba- 
bility ks (Qs) dt, the wavefuction will react and form singlet product and with probability 
ks (Qt) dt the wavefunction will be "measured" and forms a pure triplet, hence fails to react 
and escapes the singlet recombination channel as a pure triplet. Analogously, with proba- 
bility kr (Qt) dt the wavefuction will react and form triplet product and with probability 
kf (Qs) dt the wavefunction will be "measured" and form a pure singlet. Expanding (41) 
and keeping only first order terms in dt, one obtains 



p(t + dt) = p(t) - i[H, p] - tf{j+j u p(t)}dt + ^{ J 2 + J 2 , p(t)}dt + k s Jip(t)Jtdt + k T J 2 p(t)J^dt 



-!f(J 3 p(t) + p(t)J 3 - 2J 3 p(t)J+)dt - !f (J 4 p(t) + p{t)J 4 - 2J 4 p(t)J+)dt 



2 + P(t) J-3 ~ 2 J 3 P(*)4j > 

projecting into the S-T basis and making use of the relation J 3 = 1 — J 4 , equation (9) is 
recovered. 

B. Derivation of the Traditional Master Equation 

The phenomenological master equation can be derived in a way similar to the Jones-Hore 
master equation, but using different jump operators. 



\iP(t + dt)) = \,f>(t)) p=l-k s (Q s ) dt - k T (Q T ) dt 

J($\ e iH °ff dt e- iH -ff dt 

|# + 4 = |5 } = ^M P=k s (Qs)dt 

v \Qs) 

m + dt)) = \T ) = J2 j^ p=k T (Q T )dt (43) 

V \Qt) 

H ef f = H- i^Qs - i^Qt (44) 

The interpretation of (1431) is that with probability ksdt, the wavefunction attempts to 
react via the singlet channel and with probability krdt, the wavefunction attempts to react 
via the triplet channel. With probabilities kg (Qs) dt and kx (Qt) dt a reaction occurred. 
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However, nothing can be said about reactants that attempt to react but fail to do so. 



Therefore, a failed reaction doesn't return a pure singlet or triplet 

In a qualitative sense, the difference between the two master equations can be described 
using transition state theory. With a perfectly penetrative barrier, reaching the barrier but 
not reacting necessarily implies that one is in the wrong spin state, thus the fact that the 
molecule is in the singlet channel yet does not react is physically significant. With a partially 
reflective barrier, no reaction doesn't imply the reactant is in the wrong spin state - reactant 
that is in the correct spin state and reaches the barrier can still get reflected. 

Making the analogy between spin-selective chemical reaction and quantum optics, the 
difference between the quantum measurement master equation and the traditional master 
equation corresponds to measurement of the state of a quantum optical system via a flu- 
orescence detection experiment or by observing a spontaneous decay process (Fig{T]). In 



fluorescence detection 



13] . light is shone on the atom to excite it selectively from one of its 



two ground states, the bright state, into a third excited state, whereupon it spontaneously 
emits a photon and returns to the original state. The other ground state, the dark state, is 
not excited by the incident light. Hence if one observes no photon emitted after exciting the 
system, one is sure that the system is in the dark state and thus null measurement effectively 
collapses the wavefunction 



14, 



15]. 



In the spontaneous emission scenario, two excited states are coherently interconverting 
and only one state can decay to a ground state. If one observes a photon, a measurement is 
performed and the wavefunction collapsed onto the bright state before emitting a photon. 
However, not observing a photon is not physically significant and does not collapse the 
wavefunction. Of course, not observing a photon for infinitely long time implies that the 
system is in the ground state, as explained in the previous section. 
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o 




(a) Fluorescence detection experiment (b) Spontaneous emission experiment 

Figure 1: Analogy between a quantum optics experiment and a spin-selective reaction. The 
blue arrow indicates excitation to an excited state, the red arrow indicates radiative decay 
and the green arrow indicates coherent interconversion between the energy states. The white 
and black spots indicates dark and bright states respectively. 

Acknowledgments 

The author would like to thank Professor Peter Hore and Professor Jonathan Jones for 
their expert advice. This work has been supported by a Nuffield Undergraduate Research 
Bursary. 



[1] R. Haberkorn, Mole. Phys., 1976, 32, 1491-1493. 
[2] I. K. Kominis, Phys. Rev. E, 2009, 80, 056115. 
[3] I. K. Kominis, Phys. Rev. E, 2010, 81, 029901(E). 
[4] H. Carmichael, Phys. Rev. A, 1997, 56, 5065-5099. 



12 



[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 

[17 

[18 



K. Molmer, Y. Castin, J. Dalibard, J. Opt. Soc. Am. B, 1993, 10, 524-538. 

I. K. Kominis, Phys. Rev. E, 2011, 83, 056118. 

J. A. Jones, P. J. Hore, Chem. Phys. Lett, 2010, 488, 90-93. 

K. L. Ivanov, M. V. Petrova, N. N. Lukzen, K. Maeda, J. Phys. Chem. A, 2010, 114, 9447- 
9455. 

H. Wiseman, Quantum and Semicl. Opt., 1996, 8, 205-222. 
M. Plenio, P. Knight, Rev. Mod. Phys., 1998, 70, 101-144. 

H. Carmichael, An open system approach to quantum optics, Springer Verlag, 1991. 
J. A. Jones, K. Maeda, P. J. Hore, Chem. Phys. Lett, 2011, 507, 269-273. 
J. Volz, R. Gehr, G. Dubois, J. Esteve, J. Reichel, Nature, 2011, 475, 210-213. 
M. Porrati, S. Putterman, Phys. Rev. A, 1987, 36, 929-932. 

W. M. Rano, D. J. Heinzen, J. J. Bollinger, D. J. Wineland, Phys. Rev. A, 1990, 41, 2295. 
H. Carmichael, Statistical Methods in Quantum Optics, Springer Verlag, 2002. 
J. L. Skinner, D. Hsu, J. Phys. Chem., 1986, 90, 4931. 

W. C. Schieve, L. O. Horwitz, Quantum Statistical Mechanics, Cambridge University Press, 
2009. 

Appendix A: Microscopic derivation of the Jones-Hore Master Equation 



16J, the total Hamiltonian is the sum of system, environment and 



Following Carmichael 
interaction Hamiltonians: 



H = H sys + H env + Hi n t (Al) 
Moving to the interaction picture, the Liouville-von Neuman equation is 

^ = -i[H int ,p\ (A2) 

Formal integration gives 

t 

p(t) = p(0) - i j dt'[H mt (t'),p(t')} (A3) 
b 
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Substituting p(t) back into the right hand side of (A2) gives 



dp 
It 



t 



-i[H int (t),p(0)} ~ J dt'[H mt (t), [H int (t'),p(t')}} (A4) 
o 

Tracing out the reservoir degrees of freedom, and assuming that there are no initial 
system-reservoir correlations, i.e tr R ([H int (t), p(0)]) 

t 



dpsit) 
dt 



j dt'tr R {[H int (t), [H mt (t'),p(t')}}} (A5) 



where ps = trn(p) is the density matrix of the system. In the limit of weak coupling and an 
infinitely large reservoir, the motion of the system and reservoir remains factored throughout 
the evolution, hence 

p(t) = R p s (t) (A6) 
where R is the initial reservoir state. Substituting in (A5) gives 



dpsjt) 
dt 

Writing 



t 

J dt'tr R [[H^t),[HUt'),RoPs(t')}}} (A7) 



flint = X^ ? < ( A8 ) 



where Tj are operators belonging to the reservoir subspace and Sj are operators belonging 
to the system space, gives 



dp s (t) 



= ~ E / dt'tr R {[f t (t)s z (t), [f^SjW^psit')]]} (A9) 



dt 

1,3 o 

Expanding out the sum and recognizing reservoir correlation functions, we obtain 



dp s {t) 



t 



dt 

i,3 o 



- f dt'm)W)ps{t') - Sj{t')ps{t')m\ (fi(t)f j{t')) R + h. c (mo) 



where 



fi^f^O) =tr R [R r i (t)T,(t')] (All) 
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The system Hamiltonian, in the {S,T, Sq,Tq} basis, is 

/ u \ 



H. 



sys 



OJqJi J\ + U\ J<2 J2 



(A12) 



WjOO 

\ 0/ 

The environment is assumed to be a series of harmonic oscillators that absorbs the en- 
ergy dissipated from the system and two measurement devices, which are also modeled as 
harmonic oscillators 

Henv = ^2 W * a M a M + ^2 W k°fc,2 a k.2 + ^ h k h k + ^ C k C * ( Al3 ) 
k k k k 

The interaction Hamiltonian comprises two parts, a dissipative part (Hdi SS ) which is 
modeled as the Jaynes-Cumming Hamiltonian 



16| and a measurement Hamiltonian (H mea ) 
Hdiss = ^2(9kdi A Ji + g* k a ktl Jf) + ^(g^,^ + Qk^k'^) ( Al4 ) 

k k> 

Hmea = ^ h ^K&fc + O^h) + J? h ^(A'C+ + (3* k ,C k >) (A15) 
k k' 

The measurement device b measures the singlet occupancy and the device c measures the 
triplet occupancy Note that [H sys , H mea ] = and [H sys , H^ss] 7^ 0, implying that energy 
is conserved by the measurement process but not conserved by the dissipative dynamics. 



Looking at the form of the measurement Hamiltonian, at long times, t 
devices b and c will collapse to pointer states J± J\ and J% J2 respectively 



re measurement 



18J. Let 



9k a k,\e 



(A16) 



T B = Y J ulb k e- 1 ' 



(A17) 



(A18) 



(A19) 
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Writing out explicitly the operators in the form of (A8) 

f 1 = f + 1 ^ = J ie ^ ' 

f 2 = f A1 s 2 = J+e^ 

f 3 = f + 2 s 3 = J 2 e- iwi * 
f 4 = f A>2 s 4 = J+e** 1 * 

r 5 = r# + s 5 = j x 
r 6 = r c + s 6 = J 2 

At zero temperature, reservoirs A±, A 2 , B and C have the following correlation functions 



r(t)r(o; fl = o (A20) 

; 1+ (t)f + (t')) R = (A21) 
r + (t)f(t')) R = (A22) 

oo 

(T(t)T + (t')} R = Y^\kj\ 2 e-^W = J dw\k{u)\ 2 g{u)e- iuj{t '^ (A23) 
i o 

where kj is the generic system-reservoir coupling constant, and g(u>) is the density of states 
in the continuum limit. Substituting and noting that the S and T states are orthogonal and 
reservoirs Ai, A 2 , B and C are uncorrelated and thus statistically independent 



+ [J 2 + J 2/ ?s(t') - J 2 p s (t') 4]e— (f A , 2 (*)f+(0 



R 



+[J+J 2 p s (t') ~ JiJ2Mt')JtJ2\ (r c (t)r+(t')) R + h.c (A24) 



Substituting t = t — t' 
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<IPS{1 ] So dr[Jt JCpsit - r) - J.psit - T)J+]e-^ (f^XAt ~ r) 



+[J+J 2 ps(t -r)- J 2 ps(t - r)4}e-^ (f A)2 (*)f+ 2 (f - r)}^ 

Jips(t - r) - J+ Jip 5 (t - r) Jt Ji] (r B (t)r+(t - t) 



R 

+[J}J2Ps{t ~r)- J}J2Ps{t - r)J}J 2 ] (f c{t)t +(t - t)) r + h.c (A25) 

Performing the Markov approximation and replacing t — r terms in the density operator 
by t 

dp sit) 



+ [J 2 + J2Ps(t) - J2Ps{t)4\ So dre-^ (f A)2 (*)f£ 2 (t - r 



R 

+[JtJiPs(t) ~ 4JiPs(t)JtJi] So dr {?B(t)t + B {t - r) ' 



R 



+ [J+J 2 p s (t) ~ J}J 2 ps{t)J}J2] So dr \ V c{t)T+(t -r)) R + h.c. (A26) 
Ignoring imaginary frequency shifts, one can define 

t t 

J dre^ (f AA (tWUt ~ r)) R = j dr (f B (f)f+(f - r)) r = ^ (A27) 
o o 
t t 

j dre-^ (r A , 2 (t)r% 2 (t - r)) R = J dr (f c (t)f %{t - r)} = ^ (A28) 

o o 

= -y (JtJiPsit) + Ps(t) 4 Ji ~ UiPs(t)Jt) - y {4J2Ps{t) ~ Ps(t)J+J 2 - 2J 2 p s (t)4) 
_h + k L (j+j-^q + p^J+J, - 2J+J 1 ps(t)JiJi) (A29) 

Converting back to the Schrodinger picture and projecting into {S,T} subspace, (9) is 
recovered. 
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